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Abstract 

A sufficient integral criterion for a blow-up solution of the Hopf equations (the 
Euler equations with zero pressure) is found. This criterion shows that a certain 
positive integral quantity blows up in a finite time under specific initial conditions. 
Blow-up of this quantity means that solution of the Hopf equation in 3D can not 
be continued in the Sobolev space H 2 (1Z 3 ) for infinite time. 

1 Introduction 

In 1984 Beale, Kato and Maida [[I]] showed that sufficient and necessary condition for a 
smooth solution to the 3D Euler equations for ideal incompressible fluids on the time 
interval [0, to] is the finiteness of the integral, 

to 

J sup r |u;| dt < oo. (1) 



Here u is the vorticity, connected with the velocity field v by the standard relation: 

uj = curl v. 

If the integral (|TJ) is divergent, then the vorticity blows up (or collapses) in a finite time. 

Proof of this criterion is based on the local existence theorem ||. According to this 
theorem a smooth solution to the 3D incompressible Euler equations exists if the initial 
conditions v(0) belong to the Sobolev space H q (7Z 3 ) for q > 3 and the solution itself is 
from the class, 

veC([0,t o };H q )nC 1 ([0,ta};H q - 1 ), 
where the norm of H q is defined as follows, 

\\v\\i*y&)= (£ / fWfdr) ' , (2) 



with a being a multi- index. In accordance with 0, the violation of this property leads 
to divergence of the integral (|I|) and to collapse for vorticity. 

However, in order to use the criterion (|]), effectively one needs to have an explicit 
solution to the Euler equation, that is practically very difficult or impossible. Therefore, 
the main task is to find sufficient conditions which guarantee blow-up of the vorticity. 
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The first step in this direction would be to construct such criteria for models more simple 
than the incompressible Euler equations which is one of the main aims of this paper. 

We consider the simplest variant - the Euler equations for a compressible fluid without 
pressure (the hydrodynamics of dust) which are sometimes called the Hopf equations: 

^ + div pv = 0, (3) 

-1 + ( V • V)v = 0. (4) 

These equations can be successfully solved by means of the Lagrangian description. In 
terms of Lagrangian variables the equation (£|) describes the motion of free fluid particles. 
For this reason breaking is possible in this model: it happens when trajectories of fluid 
particles intersect. In terms of the mapping, describing the change of variables from 
the Eulerian description to the Lagrangian one, the process of breaking corresponds to a 
vanishing Jacobian, J, of the mapping. This process can be considered as a collapse in 
this system: the density p as well as the velocity gradient become infinite at a finite time 
in the points where the Jacobian vanishes. By virtue of this the breaking is sometimes 
called gradient catastrophe. 

In this paper we construct sufficient integral criteria for breaking within the model 
and verify what the functional conditions of the theorem U correspond to those 
for the given system. Our main conjecture is the following: breaking in the 3D case within 
the system yields a violation of Sobolev space already for q = 2 instead of q > 3 in the 
theorem [[!]. 

To complete the introduction we would like to note that the model @, (|4]) has a lot 
of astrophysical applications. According to the pioneering idea of Ya.B. Zeldovich ||, 
the formation of proto-galaxies due to breaking of stellar dust can be described by the 
system (|3|), @) (see also the review [|]). It should be noticed also that recently || the 
mechanism of breaking was used to explain collapse in incompressible fluids when, instead 
of the breaking of fluid particles, the breaking of vortex lines can happen which results in 
a blow-up of vorticity. Therefore we believe that the results presented in this paper can 
also be useful from the point of view of collapse in incompressible fluids. 

2 One-dimensional analysis 

We begin with a one-dimensional calculations where Eq. ([|) has the form 

v t + vv x = . (5) 

Consider the integrals 

I n = J v " dx 

with integer n, where integration is assumed to be from — oo to +oo. 

It is easy to show that the evolution of these integrals is defined by the relations: 

^f = -(n-l)/ n+1 , (6) 
which can be obtained by means of ([5]) and integration by parts. 
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The relations for n = 2, 3 read: 

J v 2 x dx = - J vl dx , (7) 



d 
dt 



j t Jv 3 x dx = -2jvtdx. (8) 

These differential relations allow one to find a closed differential inequality for the 
integral I 2 = J v% dx. Applying the Cauchy-Bunyakovsky inequality to ([?[) gives the 
estimation 

§<(/^*)V <») 

Substituting, 

vl dx = -- 



2 dt 2 ' 

into (|9|) we arrive at the closed differential inequality for I 2 . 

d 2 h (dl 2 \ 2 . . 

This inequality is solved by means of change of variables 

I 2 = X a > 0, 

with unknown exponent a. We will define the exponent by requiring the absence in (TO) 
of the terms ~ X a ~ 2 Xf. Hence a = — 1 and the inequality becomes, 

X tt < 0. (ff) 

Using a mechanical interpretation of X as the coordinate of a particle and taking into 
account that the particle acceleration is negative, we can immediately conclude that the 
particle can reach origin X = in a finite time if the initial particle velocity ^(0) is 
negative. For I 2 this means a blow-up: 

Hence, by elementary calculation, we can estimate the collapse time to'- 

X(0) 7 2 (0) 



to < 



Mo) | |/»(o)| 



After multiplying ( pj) by X t (< 0) and integrating the result over time, we arrive at the 
estimation from above for I 2 : 

« () s im^ry (12) 

Thus, the blow-up of I 2 takes place if the initial velocity X t (0) is negative which is 
equivalent to the initial condition, 

ht(0) > 
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ig(0) = J vldx < 0. 



If the initial distribution of Vq(x) is symmetric with respect to x, then ^(O) = 0. However, 
as follows from (0), /3 becomes negative already at t — +0 that results in a blow-up of 
the integral I<i- 

To complete this section, let us compare the estimate (0) with the exact time depen- 
dence of hit) near the breaking point. 

In order to define this dependence, we differentiate fl5|) with respect to x that results 
in the equation, 

dv x 2 ( d d d\ 

-dT = - v - \dt = Wt +v d-x)- (13) 

After integration we have 

v' (a) 

x l + v' (a)t> li4j 

where a is the initial coordinate of the fluid particle at t = 0, Vo(a) is the initial velocity, 
the prime denotes a derivative. The denominator in (|T4|), 1 + t> a (a)t, represents the 
Jacobian J of the mapping, 

x = a + v(a)t. 

The Jacobian tends to zero for the first time at t — to, defined by 

t = mm[-v~ l (a)]>0. (15) 



Near the singular point, a = a , corresponding to the minimum, ( [TBI) J can be approxi- 
mated by the expression, 

J ~ a(io — t) + j (a — a ) 2 . 

Here 



dJ 



d 2 J 

1/to, 2 7 



da 

to> a o 



> 0. 

to,a.Q 



As a result, v x takes a singularity as r = to — t — > 0: 

w (a ) 



ar + 7(0 — a ) 2 

Hence, one can see that the contribution from the singularity to the integral I2, 



(16) 



ar + 7(0 — a ) 2 
diverges as r — > and satisfies the inequality (|Xl|) . 



V ° (a ° ] ^ ~ -4, (17) 



3 Multi-dimensional breaking 

In this section we generalize the above analysis to multi-dimensions. 

In the multi-dimensional case, instead of v x it is convenient to introduce the matrix 
U with matrix elements 

10 dxi 



4 



The equations of motion for this matrix have a form analogous to ([12]): 
where 

Our aim is now to find the inequality corresponding to (|S[). Consider two scalar charac- 
teristics of the matrix U: its trace, tr U = div v and determinant, det U. From ([18]) we 
derive the following equations for these two quantities: 

^tr U = -tr (U 2 ), (19) 
at 

4 det U = -tr U ■ det U (20) 
at 



Now we introduce the positive definite integral, 



I 



J (det Uf dr, 

which, in the ID case, coincides with I2. Due to (|20|) , we have 

J tr U (det 27) 2 dr. (21) 



dl 
~dt 



This equation generalizes the equation (|7|) to the multi-dimensional case. The second 
derivative of I will be given by the expression, 



(tr Uf + tr([/ 2 



(det Uf dr. (22) 



Applying the Cauchy-Bunyakovsky inequality to the r.h.s. of fl2TiD yields 

^ < I 1/2 ■ (^J tr(U 2 ) (detf/) 2 dr) 1/2 . 



(23) 



From (|22|) we write the integral 

| (tr Uf (det Ufdv = ^4 - J tr([/ 2 )(det f/) 2 rfr. (24) 

Now, we shall estimate the integral on the r.h.s. of this equation. We shall assume that 
all eigenvalues of the matrix U are real. Such an assumption means that the matrix U is 
close to its symmetric part, S = 1/2(27 + U T ), the so-called stress tensor (here T denotes 
transpose). In this case the antisymmetric part of the matrix U, the vorticity tensor, 
Q = 1/2(27 — U T ) [|, is small compared to S. In particular, if Q = the matrix U coincides 
with S, representing the Hessian of the velocity potential $: Uij = Sij = <9 2 $ / 'dxidxj. 

2 The vorticity tensor O is connected with the vorticity u> by the relation f2jj = hcij^uik- 
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D 

Under this assumption the trace of the matrix U, tr(U 2 ) = A?, where Aj are eigen- 

i=\ 

values of the U and D dimension, becomes positive. In this case the following relation 
between traces of U and U 2 can be easily proven: 



d I f D \' 



This inequality generates the following estimate between integrals: 

J (det Ufdr J tr (U 2 )(det Ufdr > 1 (j tr U ■ (det U) 2 d^j = -I (^j 
Substitution (^4j) and (p5|) into (§3|) gives the desired differential inequality, 



(25) 



d 2 I ( 1 \ /dA 2 , . 

Its solution is sought in power form as before: I = x a . Excluding terms proportional to 
X a ~ 2 X 2 we find that a = —D and 

*-v 

For X this results in the same inequality as (|10l): X tt < 0. The criterion for attaining the 
origin X = will be also analogous: 

X t (0) < or J t (0) > . 

Almost the same form will have the estimate for the collapse time 

t < DI{0) - 1 127) 

k< iW~ ~Jm' (27) 

where (A) is a mean eigenvalue defined in accordance with 

1 x — > v 1 



( A ) = 7jE^ = ^7/ trC/(det f/) 2 dr. 



For arbitrary .D, instead of (|TT|) , the following estimation appears for /: 



A +1 (0) 

4 Comparison with exact solution 

In order to compare the estimation ( p8|) with the exact dependence of / we have to solve 
equation (|T8|) . This solution is 

f/ = f/o(a)(l + f/ (a)t)- 1 . (29) 

Here a is the initial coordinates of a fluid particle and Uq(sl) is the initial value of the matrix 
U . By introducing the projectors P^ k > of the matrix Uq(sl) (P^ 2 = P^ corresponding to 
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each of the eigenvalues Aofc(a)), this expression can be rewritten in the form of a spectral 
expansion: 

u = J2 A ? P (k) - (30) 

k = l 1 + AQkt 

The projector P^ h \ being a matrix function of a, is expressed through the eigenvectors 
for the direct (C/o(a)^ = Xoip) and conjugated (0?7 o ( a ) = <M ) spectral problems for the 
matrix £/o(a): 

if) = 4 k ur. 

where the vectors and with different n and m are mutually orthogonal: 
Hence, the determinant of the matrix U is defined by the product, 



det U 



D \ 



7 1 - 1 - 



fe =1 1 + Aofct 



From ([50) it follows also that singularity in U first time appears at t — to, defined from 
the condition || 0: 

t = min[-l/A fc(a)], (31) 

k,a 

(compare this with P7|)). 

From (p0|) , one can see that near the singular point only one term in the sum (|30|) 
survives, 

p(n) 

u ~ X x— > ( 32 ) 

r + 7 Q/3 Aa Q Aa /3 

where the projector is evaluated at the point a = a and k = n, corresponding to 
the minimum (|3TD, t = t — t, Aa = a — a , and 

a>3 da a dap 



a=ao 



is a positive definite matrix. 

The remarkable formula (^) demonstrates that i) the matrix U tends to the degenerate 
one as t — > t and ii) both parts of the matrix U in this limit, i.e. the stress tensor S and 
the vorticity tensor Q, become simultaneously infinite (compare with M). It is interesting 
to note that at near singular time the ratio between both parts is fixed and governs by 
two relations following from the definition of the projector P: 

P s = Pl + Pl P A = PsPa + PaPs 

where P$ and Pa are respectively symmetric ("potential") and antisymmetric (vortical) 
parts of the projector P. In particular, the second relation provides the collapsing solution 
for the equation for vorticity 

-— = curl [v x uj\ 
dt 1 J 

which has the same form for both compressible and incompressible cases. It is also 
interesting to note that in the sense of the criterion (JI]), the collapsing solution (|32|) 
represents the marginal solution. 
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The asymptotic solution fl32|), far from the collapsing point should be matched with a 
"regular" solution. The corresponding matching scale can be estimated as Iq ~ 7 2//3 . This 
scale Iq alone can be taken as the size of the collapsing region for (^). This remark now 
allows one to calculate the contribution from the breaking area to the integral /. 

Substituting ( |2T| ) into I and using a change of variables, from r to a, one can get the 
expression for this contribution, 

I = c t_p_ (33) 



to 



where the constant is 

C= XondetUoH-^- 

The integral is taken over the spherical volume V with coordinate center at a = a and 
size ~ Iq. Introducing a self-similar variable £ = ar -1 / 2 , one can see that the contribution 
depends significantly on the dimension D. At D = 1 this integral behaves like 

r-V2 i n 

full correspondence with fll7f). In this case, the integral over £ is convergent at large £ 
and is not sensitive to the cut-off size Iq. In the two-dimensional geometry, however, the 
integral ( |33|) has a power dependence on r, but a logarithmic dependence on l arising 
from integration on £: 

i i l ° 

J ~ lo g;^i7l 

that satisfies the inequality (|28|). 

In the three-dimensional case, the integral diverges at large scales as the first 
power of £, becomes proportional to the size of collapsing area l as r — > 0: 

J ~ Z . 

This result for = 3 formally contradicts to the blow-up sufficient condition found above. 
This contradiction indicates only that in the three-dimensional case the blow-up of the 
integral / has no universal behavior near the singular time which should be expected 
following to the universal asymptotics (|32|). 



5 Concluding remarks 

Thus, the initial condition, 

>>°- 

represents a sufficient integral criterion for the collapse if the vorticity matrix U is small 
in comparison to the stress tensor. Under this condition, the integral I = J(detU) 2 dr 
becomes infinite in a finite time. 

In turn, divergence of I means that a solution can not be continued in the correspond- 
ing functional space. In the ID case this is the Sobolev space H 1 (1Z). For D = 2, from 
the inequality, 

I< J \U\'d 2 r {\U\ 2 = Ul\ 
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together with the embedding Sobolev inequality, 

(/ \U\ A d 2 v\ < c\\u\\ mm 

it follows that ||£7||#i(-ft2) — > oo as / — > oo. In terms of the velocity, this means that the 
solution is not continued in the Sobolev space H 2 (1Z 2 ). 

In the 3D case it is possible to write the following set of inequalities, 

1/6 



/ 1/6 < (/ \U\ 6 d 3 r 



< C\\U\ 



where the second inequality represents the partial case of the Sobolev embedding inequal- 
ity 0. 



Hence, one can see that / — > oo is equivalent to the divergence of the norm (^[) for the 
Sobolev space H 2 (7Z 3 ). Thus, the requirements for strong solutions in the hydrodynamic 
model (^|), (^) are different from those for the 3D Euler equation for incompressible fluids. 
At the moment it is hard to say whether the results presented in this paper contradict 
to the theorem (see also M) or not. In any case, it is a very interesting question. It 
should be added that collapse in incompressible fluids might happen through breaking 
of vortex lines (there are some arguments both analytical and numerical |5|, || in a favor 
of such point of view). In this case the corresponding norms would also blow up for the 
same Sobolev space H 2 (1Z 3 ). 
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